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ABSTRACT 

We use a cosmological chemodynamical simulation to study how the group environment impacts 
the star formation properties of disk galaxies. The simulated group has a total mass of M ~ 8 x 
10^^ Mq and a total X-ray luminosity of Lx ~ 10^^ erg s^^. Our simulation suggests that ram 
pressure is not sufficient in this group to remove the cold disk gas from a Vrot ~ 150 km galaxy. 
However, the majority of the hot gas in the galaxy is stripped over a timescale of approximately 1 
Gyr. Since the cooling of the hot gas component provides a source for new cold gas, the stripping 
of the hot component effectively cuts off the supply of cold gas. This in turn leads to a quenching 
of star formation. The galaxy maintains the disk component after the cold gas is consumed, which 
may lead to a galaxy similar to an SO. Our self-consistent simulation suggests that this strangulation 
mechanism works even in low mass groups, providing an explanation for the lower star formation rates 
in group galaxies relative to galaxies in the field. 

Subject headings: galaxies: kinematics and dynamics — galaxies: evolution — galaxies: stellar content 
— methods: numerical 



1. INTRODUCTION 

Understanding the role environment plays in determin- 
ing the star formation properties of galaxies remains one 
of the most important issues in galaxy evolution. Since 
the work of Hubble & Humason (I93I), many studies 
have shown that clusters of galaxies contain a higher frac- 
tion of early-type galaxies (e.g., Oemler 1974; Dressier 
1980) and fewer star forming galaxies (e.g., Balogh et al. 
1997; Goto et al. 2003) than the field. Several physical 
mechanisms have been proposed to explain this trend. 
The most promising mechanisms fall into one of three 
broad categories: (1) mergers and interactions with other 
galaxies (e.g., Holmberg 1941; Toomre & Toomre 1972; 
White 1978) and/or the cluster potential (e.g., Byrd & 
Valtonen 1990; Moore et al. 1996; Bekki 1999); (2) ram- 
pressure and/or viscous stripping of the cold gas (e.g., 
Gunn & Gott 1972; Nulsen 1982; Farouki & Shapiro 
1980; Fujita & Nagashima 1999; Abadi et al. 1999; Quilis 
et al. 2000); (3) "strangulation" in which the warm and 
hot gas in the halo is stripped cutting off the supply of 
cold gas (e.g., Larson et al. 1980; Balogh et al. 2000; 
Bekki et al. 2002). 

Recent observations suggest that star formation is also 
suppressed in some groups of galaxies (e.g., Zabludoff & 
Mulchaey 1998; Wilman et al. 2005; Weinmann et al. 
2006). While galaxy interactions and mergers are often 
assumed to play a dominant role in driving galaxy evolu- 
tion in galaxy groups (e.g., Hickson et al. 1977; Carnevali 
et al. 1981; Ishizawa et al. 1983; Barnes 1989), the recent 
observations of both X-ray and HI tails in nearby groups 
(e.g., Davis et al. 1997; Sivakoff et al. 2004; Machacek 
et al. 2005; Rasmussen et al. 2006; Bureau & Carignan 
2002) suggests other mechanisms may also be important. 
Unfortunately, most studies have concentrated on rich 
galaxy clusters, and much less is known about the mech- 
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anisms important in lower mass groups (e.g., Fujita 2004; 
Roediger & Hensler 2005; Hester 2006; McCarthy et al. 
2007). In this Letter we use a self-consistent cosmological 
chemodynamical simulation to demonstrate that stran- 
gulation can effectively suppress star formation in low 
mass groups. 

2. METHOD 

The simulations were carried out using our original 
galactic chemodynamics code, GCD+ (Kawata & Gibson 
2003a,b; Kawata et al. 2006; Kawata & Rauch 2007). 
GCD+ is a three-dimensional tree A^-body/smoothed par- 
ticle hydrodynamics (SPH Lucy 1977; Gingold & Mon- 
aghan 1977) code (Barnes & Hut 1986; Hernquist & Katz 
1989; Katz et al. 1996) that incorporates self-gravity, 
hydrodynamics, radiative cooling, star formation, su- 
pernova (SN) feedback, and metal enrichment. GCD+ 
takes into account chemical enrichment by both Type II 
(Woosley & Weaver 1995) and Type la (Iwamoto et al. 
1999; Kobayashi et al. 2000) SN and mass loss from 
intermediate-mass stars (van den Hoek & Groenewegen 
1997), and follows the chemical enrichment history of 
both the stellar and gas components of the system. The 
code also includes non-equilibrium chemical reactions of 
hydrogen and helium species (H, H+, He, He+, He++, 
H2, H2", H~) and their cooling processes, following the 
method of Abel et al. (1997), Anninos et al. (1997), and 
Galh & Palla (1998). We use the UV background spec- 
trum suggested by Haardt & Madau (2001), and take 
into account radiative cooling and heating due to heavy 
elements (Raymond & Smith 1977; Cen et al. 1995). 

The system which we focus on here is one of the simu- 
lated groups of galaxies in our "Virtual Group Catalog" 
(VGRC) project. The VGRC is an on-going numerical 
simulation campaign with a goal of simulating more than 
30 groups with a range of masses and formation histo- 
ries at high enough resolution to resolve the structures of 
the member galaxies. All groups in the VGRC are identi- 
fied in cosmological simulations with a ACDM cosmology 
(Oo=0.24, Ao=0.76, f7b=0.042, h ^ 0.73, as = 0.74, and 
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Fig. 1. — Evolution of the distribution of dark matter density from an arbitrary but fixed projection [top rows), and the close-up 
rest-frame i?-band images of the target galaxy from the face-on {middle rows) and edge-on (bottom rows) views. The contours correspond 
to the cold gas density of log(p(.old/M0pc~^) = 0.2, 0.6, 1.2, 1.6. The arrows indicate the velocity of the target galaxy with respect to the 
velocity of the group. The length of the arrows of 1 kpc correspond to the speed of 20 km s~^. Bruzual & Chariot (2003) model is adopted 
for the stellar population synthesis. 
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General properties of the group at z = and the target galaxy at z = 0.51 
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^Thc dofinition of Kitayama & Suto (1996) 

^Thc X-ray luminosity in energy range between 0.01-100 keV from the extrapolation of the modeled flux from the X-ray spectrum within 
rsoo — 244 kpc, where rsoo is the radius of a sphere containing a mean density of 500 times the critical density of Pcrit — SHq/SttG at 
z = 0. 

^Thc peak gas rotation velocity around the radius of 7 kpc. 



rig = 0.95) consistent with the measured parameters from 
three-year Wilkinson Microwave Anisotropy Probe data 
(Spergel et al. 2007). We use a multi-resolution technique 
to achieve high-resolution in the regions of the identified 
groups, including the tidal forces from neighboring large- 
scale structures. The initial conditions for the simula- 
tions are constructed using the public software LINGER 
and GRAFIC2 (Bcrtschinger 2001). Gas dynamics and 
star formation are included only within the relevant high- 
resolution region, the surrounding low-resolution region 
contributes to the high-resolution region only through 



gravity. For this particular group simulation, the size of 
the high-resolution region is about ~4.5 Mpc at 2: = 0, 
and the low-resolution region is a sphere with diameter 
of about 35 Mpc. Consequently, the initial conditions 
consist of a total of 1,234,583 dark matter (DM) par- 
ticles and 783,424 gas particles. The mass and soften- 
ing lengths of individual gas (DM) particles in the high- 
resolution region are 4.41 x lO*" (2.08 x 10'') Mq and 
0.96 (1-61) kpc, respectively. The high-resolution region 
is chosen from a low-resolution simulation as the region 
within 4 times the virial radius of the identified group. 
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Fig. 2. — Top: The virial radius of the group (dotted) and the 
distance of the target galaxy from the group center (solid) as a 
function of redshift. Second: The star formation rate (black solid), 
the net mass-accretion (dotted), and mass-loss (gray solid) rates 
of baryon within 30 kpc of the target galaxy as a function of red- 
shift. Third: Relative speed of the target galaxy to the group as 
a function of redshift. Bottom: The intra-group gas density at 
the position of the target galaxy (solid) and the hot (dotted) and 
cold (gray solid) gas densities of the target galaxy as a function 
of redshift. We define Pigm the volume weighted mean den- 
sity of the hot gas between the radius of 50 and 100 kpc from the 
target galaxy. The hot (cold) gas density indicates the volume 
(mass) weighted mean density within 30 kpc. Note that the cold 
gas density is scaled by 1/100. 
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Fig. 3. — The estimated restoring force, az\d^z/dz\, for the cold 
(filled circles) and hot (open circles) gas particles as a function of 
their radius at z = 0.26, where cr^ and d^z/dz are respectively 
the surface gas density and the restoring gravitational acceleration 
along z direction which is set to the direction of the relative velocity 
of the target galaxy to the group. The horizontal line shows the 
estimated ram-pressure at the position of the target galaxy at z = 

0. 26. Note that the radius is the three-dimensional radius, and 
the restoring force is estimated at the position of each particle, 

1. e., we do not assume any symmetry, but use directly the particle 
distribution. Therefore, particles with the same radius can have 
different restoring force. The cold gas particles at large radii (r > 
10 kpc) correspond to the tail of the cold gas seen in Figure 1. 

The general properties of the group are summarized in 
Table 1. The simulation starts at z = 29.7, and the 
initial temperature and the fractions of hydrogen and 
helium species are calculated by RECFAST (Seager et al. 
1999, 2000). We turn on the UV background radiation 
at z = 6 (e.g., Fan et al. 2001). ModeUng of SN feedback 
is the most uncertain part in numerical simulations of 
galaxy formation. In the VGRC project, we have chosen 
a relatively strong feedback model, where each SN yields 
the thermal energy of 3 x 10^^ erg. 

3. RESULTS AND DISCUSSION 
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Fig. 4. — The history of the mass fraction of the hot gas (dot- 
ted), cold gas (dashed), stars (solid), and the escaped (dot dashed) 
components which are originally from the hot (upper) and cold gas 
(lower) particles within 30 kpc of the target galaxy ai z = 0.26. 
The mass fractions for the upper and lower panels are normalized 
by the mass of the hot and cold gas at z = 0.26. The particles 
which are further than 30 kpc from the center of the target galaxy 
are defined as escaped. Note that the cold gas particles can be- 
come hot gas due to heating by SN feedback. The mass of stellar 
particles decreases and the gas particle mass can increase because 
of the mass-loss from stars. 

The results of the group simulation are summarized 
in Figures 1 and 2. The top panel of Figure 1 shows 
the evolution of the dark matter density of the simulated 
group since z = 0.51. To study how the star formation 
properties of galaxies are impacted by the group, we fol- 
low the evolution of a disk galaxy on its first passage 
through the group. The properties of this galaxy (here- 
after referred to as the target galaxy) before it enters the 
group are summarized in Table 1. The middle and bot- 
tom panels in Figure 1 show face-on and edge-on views 
of the target galaxy, respectively. The contours give the 
surface density of the cold gas mass and the arrows show 
the velocity of the target galaxy with respect to the group 
velocity. Throughout the Letter, we define the cold (hot) 
gas as the gas with logT < 4.3 (logT > 4.3). 

The top panel of Figure 2 displays the time evolution 
of the distance of the target galaxy from the center of 
the group. From this panel it can be seen that the target 
galaxy first falls into the group (i.e. within the virial 
radius) around z — 0.31 and makes its closest approach 
to the group center at z ^ 0.2 before moving beyond the 
virial radius again by z = 0. The second panel in Figure 2 
shows the star formation rate (SFR) in the target galaxy. 
Particularly noteworthy is the sharp drop in the SFR 
after the galaxy passes the group pericenter, with star 
formation effectively ceasing by z ^ 0.1. The sharp drop 
in the SFR after the galaxy enters the group suggests the 
group environment can effectively quench star formation 
even in relatively low mass groups. 

After the galaxy enters the group, the cold gas dis- 
plays a clear tail morphology in the direction opposite of 
the galaxy's motion (see the z — 0.26 panel in Fig. 1). 
Such cold gas tails are often cited as evidence for ram- 
pressure stripping. To examine this possibility in our 
simulated group, we track the relative speed, v, the hot 
gas density of the intra-group medium (IGM), pioM, at 
the position of the target galaxy as a function of redshift 
(Fig. 2), from which the ram-pressure. Pram = Pigm^^j 
can be calculated. In Figure 3, we compare the estimated 
ram-pressure at z = 0.26 with the restoring force for the 
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cold and hot gas particles within the target galaxy. Al- 
though the effect of ram-pressure, which should be time 
and spatial dependent, is more complicated than in the 
simulation, the figure provides a qualitative comparison 
and supports our results. Figure 3 demonstrates that the 
ram-pressure experienced by the target galaxy is well- 
below the amount required to strip the cold gas. We also 
examine how the fraction of the hot gas and cold gas par- 
ticles within 30 kpc (an arbitrary choice to include the 
majority of the cold gas but avoid contamination from 
the IGM) of the target galaxy evolves from z = 0.26 to 
z = 0.0. Figure 4 shows that virtually none of the cold 
gas escapes. Interestingly, even the particles in the cold 
gas tail at z = 0.26 do not escape from the galaxy. Thus, 
we conclude that ram-pressure stripping of the cold gas 
component is not the cause of the downturn in the SFR. 

However, the majority 90 %) of the hot gas is indeed 
stripped since z = 0.26. In addition, a significant fraction 
of the hot gas is converted into cold gas. This suggests 
that hot gas halos are potentially an important source 
of cold gas in disk galaxies. Furthermore, since ram- 
pressure is strong enough to strip the hot gas (Fig. 3), 
the new supply of cold gas is cut-off. By z ~ 0.1 almost 
all the cold gas has turned into stars and star formation 
ceases. This is also demonstrated in the second panel 
of Figure 2 which shows the net baryon mass-accretion 
and mass-loss rates of the target galaxy, i.e., simply the 
difference in the total baryon mass divided by the time 
between the simulation outputs. Before the target galaxy 
falls into the group, the mass-accretion rate is compara- 
ble to the SFR. However, the mass-loss due to the hot 
gas stripping becomes significant around z ^ 0.26, and 
the mass-accretion rate is no longer sufficient to sustain 
the SFR. Therefore, we conclude that star formation is 
quenched in the target galaxy because of "strangulation" 
(Larson et al. 1980). 

Figure 2 shows some enhancement of star formation 
since z ~ 0.26. This suggests that the group environment 
somehow enhances star formation, which accelerates the 
strangulation process. The bottom panel of Figure 2 
displays the hot and cold gas density in the target galaxy. 
The hot gas density decreases since z ^ 0.26 due to the 
stripping and becomes close to the IGM density. On the 
other hand, the density of the cold gas where stars form 
significantly increases. One possibility is that the cold 
gas is compressed due to ram-pressure, which leads to 
enhanced star formation (see also Schulz & Struck 2001; 
Bekki & Couch 2003). 

Recently, Agertz et al. (2007) have shown that the SPH 
scheme can underestimate fluid instabilities, such as the 
Kelvin- Helmholtz instability (KHI), which may be im- 
portant for the stripping of cold gas (Nulsen 1982; Mori 
& Burkert 2000; Mayer et al. 2007). To test the impact 
of the KHI in our simulation, we follow Mori & Burk- 
ert (2000) to estimate the timescale of the KHI (i.e., 
TKH = Mgas(< r) / {nr'^ pigmv)) for the target galaxy at 
z = 0.26 when the hot gas stripping becomes signifi- 
cant. At r = 10 kpc (an arbitrary choice close to a 
radius where the hot gas stripping is expected in Fig. 3), 
we find that tkh ~ 2.8 Gyr (Mgas(< r) = 1.6 x 10^ 
M0, piGM = 2.84 X 10~^* g cm-^, and v = 416 km s~^ 
at = 0.26). This is longer than the timescale of the 
strangulation process for the target galaxy, i.e., between 



z = 0.26 and 0.1 (- 1.7 Gyr). Hence, the KHI is not 
the dominant mechanism of the gas stripping for the tar- 
get galaxy (see also Mayer et al. 2006; McCarthy et al. 
2007). Note that the effect of KHI is more complicated 
than our simple estimate. Therefore, our simulation may 
underestimate the stripping mass. 

It is worth noting that at z = the target galaxy is 
once again beyond the virial radius of the group, and 
still has no star formation. Although the target galaxy 
has consumed its cold gas disk via star formation by 
z = 0, the disk structure of the stellar component is 
still intact (Fig. 1). Therefore, the properties of the tar- 
get galaxy are morphologically similar to those of an SO 
galaxy. Thus, we suspect that strangulation may pro- 
duce a populations of SO galaxies (Larson et al. 1980) in 
and around poor galaxy groups. The current simulation 
suggests it would be interesting to more quantitatively 
compare a statistical sample of simulated galaxies with 
observed SO galaxies in the future. 

Recent observations of the spiral galaxy NGC 2276 in 
the NGC 2300 group support our conclusion that stran- 
gulation is important in galaxy groups. Rasmussen et al. 
(2006) have performed a detailed study of this system 
with Chandra and conclude that ram-pressure stripping 
is not sufficient to strip the cold gas in the disk. However, 
NGC 2276 is losing gas from the hot phase at a rate of 
^ 5 Mq yr^"'^. Given this mass loss rate, the galaxy will 
likely exhaust its current supply of cold gas within 1-2 
Gyr. Similar to what we find for our target galaxy, Ras- 
mussen et al. (2006) conclude that once star formation 
has ceased, NGC 2276 will have properties similar to an 
SO galaxy. The NGC 2300 group is more massive than 
the group in our present study. Our results suggest that 
evidence for strangulation may be found in lower mass 
systems. 

The simulated group is a fairly relaxed system with a 
dominant central elliptical galaxy. Strangulation is likely 
to be most effective in dynamically evolved groups be- 
cause of their enhanced IGM. Furthermore, galaxies en- 
tering relaxed groups will be subjected to stripping over 
larger path lengths than galaxies entering dynamically 
young systems which consist of individual halos and not 
a common envelope. The presence of a significant X- 
ray emitting IGM is a good indicator that a group is 
relaxed (Mulchaey 2000). As nearly all X-ray groups 
contain a dominant central elliptical galaxy, while sys- 
tems dominated by spirals tend to be non-detected in 
X-rays (Mulchaey et al. 2003; Osmond & Ponman 2004), 
strangulation is likely more effective in systems with a 
central dominant elliptical galaxy. If true, this might 
provide an explanation for the recent finding by Wein- 
mann et al. (2006) that the early-type fraction is higher 
in systems with an early-type central galaxy than in halos 
with a late- type central galaxy ("galactic conformity"). 
The larger sample of the simulated groups in the VGRC 
will allow us to test this idea. 

Although this Letter highlights only one case of stran- 
gulation, we are seeing a similar suppression of star for- 
mation in other groups in the VGRC. With the large 
range of properties among the VGRC, we hope to be 
able to quantify the importance of strangulation and de- 
termine how the suppression of star formation depends 
on factors such as the galaxy's mass and orbit. The sim- 
ulated groups in the VGRC will also allow us to explore 
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the importance of other mechanisms such as galaxy in- 
teraction in driving galaxy evolution in groups. 
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